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O . ABSTRACT 

o ■ 
<N : 

pH . Beginning with the 7-ray bursts GRB 970228 and GRB 970508, a standard model 

, ^ , for the interpretation of GRB afterglows emerged involving synchrotron emission from 

, a constant energy blast wave expanding into a constant density, "interstellar" medium. 

CSJ I However, a massive star origin for GRBs implies a stellar wind environment, probably 

' a Wolf-Rayet star, and we have previously suggested wind interaction models for 

> ■ the afterglows of GRBs 980326, 980519, and 980425/SN 1998bw. Here, we extend 

■ the theory of afterglows in winds, considering strong cooling phases, the transition 

■ to nonrelativistic motion, jets, and prompt, reverse shock emission. We find that, 

, compared to the interstellar case, the optical prompt emission in the wind case could 

Q\ I have a comparable magnitude, but would die off faster. We examine the afterglows 

I of other well-observed GRBs in the context of wind interaction models, and find that 

i-^^l GRBs 970228 and 970508 are likely wind interactors. The revision in the nonthermal 

I ' afterglow emission from GRB 970228 caused by the recognition of late supernova 

2 ' emission favors wind interaction. The radio evolution of GRB 970508 provides 

c/3 . especially strong evidence for wind interaction. For GRB 970508, the observations 

• • , suggest a density that is compatible with that expected in a Wolf-Rayet star wind. 
^ , 

• w^ , Finally, observations of the afterglow evolution of GRBs 990123 and 990510 and the 

^ I prompt optical emission from GRB 990123 favor interstellar interaction models, which 

' would suggest compact star merger progenitors for these objects. 



Subject headings: gamma-rays: bursts — stars: mass loss — stars: supernovae: general 



1. INTRODUCTION 

The fireball model for GRBs (gamma-ray bursts) led to predictions of the afterglow emission 
that might be expected when the energetic shock wave encountered the surrounding medium 
(Katz 1994; Meszaros, &: Rees 1997). The subsequent optical and X-ray observations of the 
afterglow from GRB 970228 appeared to confirm the predictions of the simplest afterglow model 
(Waxman 1997a; Wijers, Rees, & Meszaros 1997). This model involved synchrotron emission from 
electrons accelerated to a power law energy spectrum in a relativistic blast wave expanding into 
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a constant density, presumably interstellar medium (ISM). In particular, the expected relation 

between the flux spectral index and the power law rate of flux decay was in approximate accord 
with the observations. This model has become the "standard model" for the interpretation of 
GRB afterglow observations. It was used to make predictions of bright optical emission in the 
early phases when a reverse shock front is present (Sari k. Piran 1999b). The observation of a 
bright flash from GRB 990123 (Akerlof et al. 1999) gave basic confirmation of this aspect of the 
model (Sari &: Piran 1999a; Meszaros & Rees 1999). The expectation of jets for the initial energy 
deposition led to predictions of the effects on the light curve as the jet slowed (Rhoads 1997, 1999; 
Sari, Piran, &: Halpern 1999). The observations of the afterglow of GRB 990510 confirmed the 
basic features expected for jet deceleration (Harrison et al. 1999). 

These successes of the standard model give confidence that it is essentially correct. However, 
there has been increasing evidence that at least some GRBs have massive star progenitors, 
as initially suggested by Woosley (1993) on theoretical grounds. Paczyhski (1998) noted that 
the available evidence on the location of GRBs in their host galaxies indicated a link to star 
formation. A more direct link to massive stars is provided by the presence of a supernova, the 
light from the exploded star matter. GRB 980425 was probably associated with the relatively 
nearby Type Ic supernova SN 1998bw (Galama et al. 1998d; Kulkarni et al. 1998), and evidence 
for supernova type emission has now been found in GRB 980326 (Bloom et al. 1999b) and GRB 
970228 (Reichart 1999; Galama et al. 1999b). 

The importance of a massive star origin for the afterglow evolution is that the GRB blast wave 
should be expanding into the stellar wind of the progenitor star. Dai k. Lu (1998), Meszaros, Rees, 
& Wijers (1998) and Panaitescu, Meszaros, & Rees (1998) described some features of afterglow 
evolution in a p oc stellar wind. Chevalier & Li (1999, hereafter CL) made specific estimates 
for expansion into the wind of a Wolf-Rayet star. Li &; Chevalier (1999) and CL identified GRB 
980425/SN 1998bw, GRB 980326, and GRB 980519 as likely circumstellar wind interactors based 
on their afterglow evolution. In the cases of GRB 980425 and 980326, this association is supported 
by the presence of supernova emission. In GRB 980519, a supernova would have to be somewhat 
fainter than SN 1998bw (Bloom 1999). CL identified GRB 990123 as a likely interstellar interactor 
based on its afterglow evolution. GRBs of this type are not expected to be accompanied by 
supernovae. 

CL placed the afterglows of GRB 980519 and GRB 980326 in the wind category partly based 
on their relatively rapid rates of decline of optical emission. An alternative explanation for the 
decline is that the emission is from a laterally expanding jet (Sari, Piran, &: Halpern 1999; Halpern 
et al. 1999). However, in the case of GRB 980519, CL were able to use radio observations made 
during the first 3 days (Frail et al. 1998b) to further constrain the model and to predict the radio 
evolution. Radio data on GRB 980519 over the first 63 days are now available (Frail et al. 1999) 
and they are in agreement with the prediction of CL. The jet model does not appear to fit the 
data as well, but scintillation effects cause sufficiently large uncertainties in the radio fiuxes that 
the jet model is also acceptable. 
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In view of the increasing evidence for massive star progenitors of GRBs, our aim here is 
to further develop the theory of wind interaction and to examine the data on afterglows in the 
context of wind and interstellar interaction models. In § 2, we extend the theory of afterglows in 
winds to consider strong cooling, nonrelativistic evolution, and jets. Our discussion is guided by 
previous discussions of the constant density interaction case, and we contrast the two situations. 
The differences are expected to be especially significant for prompt emission because of the large 
difference in the ambient density at early times. We treat prompt, reverse shock emission for the 
wind case in § 3. In § 4, we examine the data on the best observed afterglows in the context 
of the wind and constant density models. We find that wind models can explain a number of 
observations previously regarded as puzzling, pointing to wind interaction models for these cases. 
Although wind models are indicated for some observed afterglows, there are others that are better 
described by constant density interaction. We discuss some implications of this result in § 5. Our 
conclusions are listed in S 6. 



2. AFTERGLOW LIGHT CURVES IN WINDS 

2.1. Blast Wave Hydrodynamics 

The basic model for GRB afterglow hydrodynamics involves a relativistic blast wave expanding 
into the surrounding medium (e.g., Meszaros k, Rees 1997). During the very early evolution, the 
GRB ejecta play a role and we treat the hydrodynamics of that phase in § 3. The later blast 
wave, which is dominated by the energy deposited in the external medium, can be described by a 
self-similar solution (Blandford &: McKcc 1976). The solution for expansion in a constant density 
medium has been widely used in GRB studies and here we discuss some of the basic results for a 
medium with density p = Ar~'^ , where A is a constant. The main problem is to determine the 
blast wave characteristics seen by an external observer. 

For an ultrarelativistic, adiabatic blast wave, Blandford & McKee (1976, their equation [69]) 
find that 

^= 17-4S ' 

where E is the total energy, F is the Lorentz factor of the shock front, and R is the shock wave 
radius. The ultrarelativistic shock condition yields the Lorentz factor of the gas, 7 = F/v^. 
Because E is constant, we have the standard result 7 oc R~^^~^^/^. For an observer viewing the 
blast wave along the line of sight to the center (not at a cosmological distance), consideration of 
emission from the blast wave at two times in the blast wave frame yields the time in the observer's 
frame (e.g.. Sari 1997; Panaitescu & Meszaros 1998; Dai & Lu 1998): 



Rl 



4(4-s)7ic' 

where the subscript L refers to the line of sight. For s = 0, we have t = RL/{16jj^c), as in Sari 



(2) 
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(1997). For the wind case considered here (s = 2), we have t = Ri/iS'jlc) (Dai k Lu 1998). The 
blast wave is undecelerated for s = 3 and we obtain t = Rl/ {A'^\c)] this can also be expressed as 
t = Ri/(2T'j^c), which is the standard result for a point moving toward the observer at constant 
velocity. Substitution into equation (|^) then yields 



Rl 



(4-s)(17-4s)i?t]^/(^-^) 



AttAc 



We then have Rl = {17 Et / t: pocY ^ ^ ■, where po is the ambient density, for s = and 
Rl = (9£;t/27r^c)i/2 for s = 2. 

The problem with these expressions is that they apply only to the line of sight. In observing 
a burst, the emission can be dominated by emission away from the center because the burst 
is observed at an earlier time when it may have been brighter (Waxman 1997c, Panaitescu &: 
Meszaros 1998; Sari 1998). The appearance depends on the evolution of the burst, which itself 
depends on the observing frequency. Panaitescu & Meszaros (1998) considered these issues for 
wind and constant density surroundings and for radiative and non-radiative blast waves. The 
radius and Lorentz factor of the typical material that is observed can be written as i? = C^L and 
7 = C~^''^7L) where C is a constant and the relation between R and 7 follows the 7 oc i?^^/^ relation 
expected for an adiabatic blast wave in a wind. Panaitescu &l Meszaros (1998) find C, = 0.56 at 
high frequencies where oc zy~(P~^)/^ and C, = 0.78 at low frequencies where Fy cx v^^^ . We thus 
have for a distant observer 



i?=l.lxlO^^(i±^)" Eli^A-^'h^cu. (4) 

and 



^=^K^)'"^-*-^*"*'i" (5) 

at high frequencies, where E^2 is the blast wave energy in units of 10^^ ergs, tday is the observer's 
time in units of days, A = M^/AttV^ = 5 x 10-*^^^^ g cm"-"^, is the mass loss rate, and 
is the wind velocity. The reference value of A corresponds to = 1 x 10~^ Mq yr~^ and 
Vw = 1000 km s~^ (see CL for a justification of these values in terms of the wind from a Wolf-Rayet 
star). In addition, the cosmological redshift z enters because of time dilation. At low frequencies, 
the coefficient in equation (Q) is increased to 1.6 x 10^^ cm and that in equation (^) is reduced to 
5.0. 

At early times, the adiabatic assumption is expected to break down, because all of the 
electron energy can be radiated by synchrotron radiation close to the shock front. The importance 
for the hydrodynamics is determined by the fraction of the total shock energy that goes into the 
electrons, eg. The blast wave is radiative either if eg is close to 1 or eg is small, but the magnetic 
and nucleon energies efficiently couple to the electron energy. In the limiting case, all of the shock 
power is radiated and we have a radiative blast wave (Blandford Sz McKee 1976; Vietri 1997; Sari 
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et al. 1998). For s = 2, the blast wave radius along the line of sight is (Bottcher & Dermer 1999) 

B -f^^V' (6) 

where Eq is the initial energy and To is the initial Lorentz factor of the shock front. The evolution 
of the shock wave Lorentz factor is given by 



(7) 



The estimated values of C for this case are 0.83 (low frequency) and 0.68 (high frequency) 
(Panaitescu &; Meszaros 1998), where now R = CRrad,L and T = C~^^rad,L- 

The limiting radiative case is unlikely to be achieved and the actual case is probably 
intermediate between the adiabatic and strongly radiative cases. If e = ee^rad is a constant, where 
the synchrotron radiative efficiency trad is between and 1, we have (Bottcher k, Dermer 1999) 

During the radiative phase, e^ad ~ 1 so that these expressions apply to the case where is 
constant. 



2.2. Afterglow Properties with Slow Cooling 

Our treatment of afterglow light curves in winds follows the discussions of afterglows in a 
constant density interstellar medium (Rees k, Meszaros 1997; Waxman 1997a,b; Sari, Piran, & 
Narayan 1998; Wijers k, Galama 1999). We presume that electrons are accelerated in the blast 
wave shock wave to a power energy distribution, TV (7) oc 7"^ for 7 > 7^ where 7^ is the minimum 
Lorentz factor at the shock front. The value of 7^ is 



where X is the hydrogen mass fraction, nip and rrif, are the proton and electron mass respectively, 
and eg is again the ratio of the energy density in electrons to the total postshock energy density. 
Our expression agrees with that of Sari et al. (1998), but differs slightly from that of Wijers k 
Galama (1999) who identify eg as the the ratio of the energy density in electrons to the postshock 
nucleon energy density. We take X = because we are assuming the wind is from a Wolf-Rayet 
star. The electron energy for which synchrotron losses are important is estimated as that at which 
the cooling time for electrons with pitch angle 7r/2 equals the expansion time: 

evrmec 

C^T7-t>t0ttl0C 
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where a^. is the Thomson cross section, Stot is the total field strength, and tioc = t/{l + z) is the 
time observed by an observer cosmologically local to the burst. 

The synchrotron spectrum of the afterglow can be divided into 4 power law sections that 
are separated at 3 characteristic frequencies: the synchrotron self-absorption frequency va^ the 
characteristic frequency Um emitted by electrons with Lorentz factor 7m,, and the frequency 
at which synchrotron losses become important Uc. (Sari et al. 1998). If the frequencies are 
ordered va < < i^c, the four sections of the spectrum can be described by F^, oc with 
P = 2,1/3, —{p — l)/2, —p/2 going from low to high frequency. The observations of GRB 970508 
on day 12 can be approximately represented by this spectrum (Galama et al. 1998b) and we 
initially assume this ordering of the characteristic frequencies in this section. The slow cooling 
condition can be expressed as < fc (Sari et al. 1998). The peak flux, Fj^,max) occurs at z/^ for 
this case. 

In any power law segment of the spectrum, the fiux evolution is expected to follow a power law 
F,^ oc v^t"'. For a constant density medium (s = 0), the standard result is a = —3{p— l)/4 = 3/3/2 
for Vm < < Vc and a = — (3p — 2)/4 = (3/3 + l)/2 for Vc < y- In a wind (s = 2), we have 
a = —{3p — l)/4 = (3/3 — l)/2 for < ^ < i^c and the same evolution as the s = case for 
Vc < V. A signature of wind interaction afterglows in the non-cooling, high frequency phase of 
evolution is a relatively rapid rate of decline. The plausible assumption that most of the electron 
energy is near 7^ requires p > 2 and a < —1.25. Many optical afterglows are observed to initially 
have — 1.3 and so are unlikely to be in the non-cooling, wind category. However, if the optical 
emission is in the cooling regime, a = — (3p — 2)/4 as in the s = case and the requirement p > 2 
implies a < —1.0. More candidate afterglows are potentially in this category, but there is the 
additional requirement of a moderately steep observed spectral index: /3 = —p/2 so that p > 2 
implies f3 < —1.0. 

The observed values of -Fv.max, i^m, I'A, and Uc have been used to find the blast wave energy E, 
ambient density n, electron energy fraction eg, and magnetic energy fraction in the context of 
s = models (Wijers &: Galama 1999; Granot et al. 1999). In wind models, the constant density 
is replaced by the wind density p = Ar~'^, where A = M^/AirV-vj = 5 x 10-*^-^^* g cm~-^, as discussed 
below equation (^). In CL, Gxprcssions W6r6 derived for the characteristic frecjuencies a-nd -^t/^max 
in terms of the model parameters. We repeat the expressions for completeness: 

F„,„„ . 20 (^^) " (l±i) (g) Eli'A.t^^ (U) 

— «'1^)"(i&)'XS)"<C^-. (-) 

... 2 X 10- (i±i) (^) E^A:',^ H., (14) 
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where our expressions assume a flat universe with Hubble constant Ho = 65 km s~^ Mpc~^. We 
also assume that the composition of the wind gas is hydrogen depleted and that the electron 
spectral index is p ~ 2.5; the variation due to different values oi p is less than other uncertainties. 
Our expression for u^. uses the estimate of Sari et al. (1998). The value of Vc is sensitive to the 
value of C discussed in § 2.1 (oc C,^)] we have taken the high frequency value. The expression 
for the flux here and elsewhere in the paper can be generalized to any cosmology by replacing 
(Vl + z — 1) by (iL/(9.23 Gpc \/l + z), where is the luminosity distance. In the case of a fiat 
universe, 

d, = ^{l + z-Vl + ~z). (15) 

We now invert these expressions in order to solve for the model parameters. The observations 
are taken to all refer to the same day, t. Then 

i^. ~ d X lU yx y^^^j y j (^^^9 y^^^^ ^^j y^^^^ ^^j ergs, 

(16) 

^= » - i^r {"mr i^r i^r^' i^r^ 
— lov- i^y' {^r {^r^ 

where 



l + z , \/2x - 1 

X = — - — and y = — ^ . 

2 ^ V2- 1 

The afterglow light curve depends on how the break frequencies and the peak flux, -^//^max? 
evolve with time and here we consider the light curve without cooling. As in Sari et al. (1998) 
for the ISM case, we can define critical times at which the break frequencies pass through a fixed 
frequency v. The time that crosses the observed frequency is 

*- - ^» (in)"' (n) ^'^'^ t^") 

where uw is u in units of lO^'^ Hz. The time that va crosses the observed frequency is 

U = 50 (l±i) (i^) {is.) E-"^.-" day. (21, 

Provided tA < tm, the light curve is as follows: F,y (x (x t ioi t < tA, (x -Fi/,maxt'm^''^ oc for 
tA<t< tm, and Fi, oc F^^^naxi^m'^^^'^ OC for tm < t. The condition that tA = tm leads to 

a critical time and frequency: 

l + z\ f eeV'/' fesV/' ^ .-4/3 



tAm = 80(^](^] (^) £52Ar''days, (22) 
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and 

-Am = 7 X 10' l^^y'l^-^y'E^.'Al Hz. (23) 
For v > UAm, the light curve is as described above. 

For u < VAmi there is a regime where the spectrum is characterized by Um < -A- In general, 
the flux in the self-absorbed part of the spectrum can be described by oc R^v'^{'^ehm), where 
7e is the Lorentz factor of the electrons that are responsible for the emission at frequency v. Then, 
for V < I'm, 1e = 7m and Fi, oc tz^^ as before. For i^m < - < t^A, {je/jm) = ('^/^'m)^''^ so that 
-Fjy oc iJl^v^l'^ . Above va, we again have Fy oc the v^l"^ part of the spectrum 

is no longer present. The evolution of v^. is always oc but now Fy^ oc t~'^ . The peak of the 

spectrum is at va oc t-[3{2+p)]/[2{4+p)] ^ and F^^ oc t-(i+4p)/[2(4+p)] , xhe light curve is thus described 
by Fy oc tv^ up to t = t^, followed by F^ oc t'l^v^l'^ and Fy oc t-{3p-i)/4;^-{p-i)/2_ ^an be seen 
from the estimated value of VAm that these considerations are relevant to radio observations of 
afterglows. 



2.3. Fast Cooling Case 

The description of the evolution given in CL and in §2.2 assumes that < J^c- Although 
synchrotron cooling is important for the high energy electrons, it is not important for the electrons 
emitting near which have most of the energy. Before some time, to, when = t^c, cooling 
of these lower energy electrons is rapid compared to the age and this phase can be referred to 
as having fast cooling (Sari et al. 1998). If eg is close to 1, the blast wave steadily loses energy 
and the hydrodynamic evolution is termed radiative. If ee is small and the electron energy does 
not couple efficiently to the ion or magnetic energy, the hydrodynamic evolution is approximately 
adiabatic even for t < to- We consider this second case to be the most likely and concentrate on 
it. For Um > I'c, tfie expressions for Fj^^max, J^m, and Uc (equations |jll|], and 0|) remain 
unchanged; the maximum flux now occurs at Vc instead of Vm- The spectrum can be described by 
Fy oc with P = 1/3, —1/2, —p/2 going from low to high frequency and the breaks at Uc and Um, 
respectively (Sari et al. 1998). However, equation (|l2|) for va is no longer applicable. 



Equations (|13|) and (|1) imply that 



to = 2[^)(^)(^]A^d.ys. (24) 



0.1/ vo.i 

For standard parameters, the transition to slow cooling occurs at a later time for the wind case 
compared to the ISM case because of the higher densities that the shock front encounters at early 
times. 



The afterglow light curve again depends on how the break frequencies and the peak flux, 
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Fu,max, evolve with time. The cooHng time for the wind case is 

= 2 X 10-5 (^)' {^)'^^2A>'o days. (25) 

The effects of coohng at a frequency i' are important for t < and fast coohng can affect the 
blast wave evolution up to time toi so that cooling has effects up to t = max(tc,to)- The time 
that Vm crosses v is given by equation (^). There are two possible orderings for the three times: 
tm < to < tc and tc < to < tm- These cases are divided by a critical frequency, Vq = i^cito) = i^mito)'- 

„„.3.10-(l±i)-'(i^)"'(iiL)-^i'^.-/^H. (.0) 

When u > h'o, the ordering tm <to < tc applies and the evolution can be called the high-frequency 
light curve (cf. Sari et al. 1998 for the ISM case). The evolution is described by oc t~^/'^v~^/'^ 

{t < tm), f-(3p-2)/4^-p/2 (tm<t< to), t-(3p-2)/4^-p/2 (^^ < ^ < Q, and t-(3p-l)/4j,-(p-l)/2 

{tc < t). When u < Uo, tc < to < tm applies and we have the low- frequency light curve: F^, oc 

^-1/4^-1/2 < r2/3i.l/3 {tc<t< to), t%^'^ {to<t< tm), and i-(3P-l)/4i.-(P-l)/2 (^^ < t). 

These light curves are distinct from the ISM case, especially because of the early importance of 
cooling. 

Self-absorption is important throughout the cooling regime {t < to) if tA > to- From equations 
(|2^ (which assumes the low-frequency light curve) and ([2^) , the frequency at which tA = to is 
given by 

= s X 10- (i±i)-' {^y'" {^y"'E;.rAi'^ h. (27) 

With typical parameters, i'ao < t^o so that the use of equation ( p2[ ) is justified. At radio 
wavelengths, the synchrotron emission is self-absorbed during the fast cooling period. At optical 
and X-ray wavelengths, self-absorption is important only during the early phases of the fast 
cooling period. 

The effects of self-absorption during the fast cooling period can be found in the same way 
as described above. The frequencies Vm and Vc are again given by the standard expressions 
(equations and |l^), and the behavior of ua can be found by considering the complete spectral 
evolution. For the cases ua < i^c < I'm and Uc < va < I'm, the high frequency properties are the 
same. For v > Vm, F^ oc t~{3p-2)/4j^-p/2 ^ below Um, we have Fy oc t~^/'^v~^/'^ . We also 

have Fp^ oc t^/^. In the lowest frequency range, where self-absorption is important, the electrons 
responsible for the emission have a Lorentz factor corresponding to peak emission at Vc- The 
evolution is thus given by Fp oc B?{uc/i'm)'^^'^ t^v"^. For ua < ^c, it is now straightforward to 
show that UA oc t~^/^, F^^ oc t~^/^ , Fi,^ oc t^^/^, and F^ oc t^'^/^u^/^ for va < v < Uc- For Vc < fA, 
we have ua oc t"^/^, F^^ oc t^, F,^^ oc t^/-^^, and Fj, oc tJl'^v^l'^ for Vc<v < va- 
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2.4. Light Curves 

The consideration of complete light curves for relativistic, spherical expansion involves the 
combination of both cooling and adiabatic evolution. At a particular frequency, cooling effects are 
important up to a time t = max(to, tc). The relevant transition times for a light curve are tA, to, tc, 
and tm- These four transition times allow many possible light curves, but our previous discussion 
allows us to limit the possibilities, on the assumption that the blast wave conditions are not far 



from the typical conditions that we have chosen. In particular, equations (22) and (|23| ) show that 
we expect tA = tm at a relatively late time and low frequency. We thus expect tA < tm at higher 
frequencies, including all those in the high frequency case with regard to radiative cooling. 

The orderings of the transition times of interest are thus: A) tA < t^. < to < t^.; B) 
tA < tc < to < tm] C) tc < tA < to < tm; D) tc < to < tA < tm] and E) tc < to < tm < tA, where 
the listing is from high frequency to low frequency light curves. The frequency dividing light curve 
A from B is Vq (equation pS]), that dividing C from D is i'Ao (equation [27|), and that dividing D 



from E is UAm (equation IQ). The frequency dividing curve B from C is that at which tc = tA- 
By referring to the previous sections, the light curves for each case can be constructed. They are: 

A) OC t7/4^5/2 < i^)^ t-1/4^-1/2 ^tA,tm), t-(3p-2)/4^-p/2 (i^^^^)^ ^-(3p-2)/4^-p/2 (^^^^^)^ 
^-(3p-l)/4^-(p-l)/2 (^^ < B) ^ ^7/4^5/2 < ^-1/4^,-1/2 ^tA,tc), t'^/^U^^ {tc,to), t^U^/^ 

{to,tm), t-(^P-^y^U-(P-^y^ {tm < t); C) OC t^^V^ (t < Q, t^U^ {tcM), t-^I^V^I^ {tA,to), 
(^t^^tm), t-(3p-l)/4^-(p-l)/2 (^^ < i). ^ ^7/4^5/2 < ^2^2 ^^2 

(.tA,tm), t-(3p-l)/4^-(p-l)/2 < i). E) F, OC t^^U^/^ {t < tc), t^V^ {tc,to), tv^ {to,tm), 

^7/4^5/2 [tm,tA), t~('^P~^)/^zy~(P~^)/2 (tA < t). The times in parentheses give the time range for 
that particular section of the light curve. The light curves are illustrated in Fig. 1. In the figure, 
tA is omitted from curve A and tc from curves D and E because they occur at very early times. 

For our typical parameters, both X-ray and optical light curves are of type A. It can be seen 



from equation ( P0|) that tm is typically in the lO's of seconds range at these high frequencies. 
The self-absorbed part of the light curve occurs considerably earlier, when the shock evolution is 
still affected by the energy deposition and our model is not applicable. The radio light curve is 
typically of type D. The light curve properties are like those described by CL, except we have now 
found that a steeper rise is expected for t < to- 

For the above light curves, we have assumed that the hydrodynamic evolution is described 
by an adiabatic blast wave. For t < to, the energy loss by electrons can affect the evolution. In 
§ 2.1, we gave the expected evolution for a fully radiative blast wave. In this case for a p oc 
ambient medium, standard expressions (Sari et al. 1998; Wijers &: Galama 1999) yield Vc t^^^, 
Vm OC t~^/^, and F;y,max OC We thus have Fi, = F,^,max(i^/i^c)"'^^^ oc t^'^/^v^/^ for ua <v < Vc, 

Fu = F,^^^^{v/vc)-^'^ OC for vc<y< ym, and F, = F,^{v /vm)~^/^ oc i-(5p-2)/6^-p/2 

for Um < V- The evolution at optical and X-ray wavelengths is most likely to be in the 
Vm < V regime; for -p = (2.0,2.5,3.0), the evolution of a radiative blast wave is given by 
Fy OC ^(-1-33,-1.75,-2.17)^ compared to F^ oc ^(-i-OO, -1.38, -1.75) foj. adiabatic blast wave. There 
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ls a significant steepening of the light curve for the radiative case. The radio Hght curve is likely 
to be in the self-absorbed regime. For v < va, J^c, we have F^, oc E?v'^{uc/vmY^'^ oc i^l'^v^ . 

The actual situation is likely to be intermediate between adiabatic and radiative evolution 
of the blast wave and we use equations to find the time dependences in that case. We find 
i/c oc t(2-^)/(4-e)^ ^ ^-(6-e)/(4-e)^ F^^rnax OC The resulting flux evolution for 

u>u^isF,oc i-[(6-.)p-2{2-e)]/[2(4-e)]^-p/2; ^hus, for 6 = 0.1, we have F, oc ^(-1-03,-1.40,-1.78) 
for p = (2.0,2.5,3.0), and for e = 1/3, we have F^ oc t(-i-09,-i.48,-i.86)^ ^g^^ ^^^^^ ^^^^^^ ^^le 
change from the adiabatic case is small for e = 0.1, but that it becomes significant for e = 1/3. 
The radio evolution is likely to be characterized hy u < v^, which is below Vc and Um. We then 
have oct(s-3^)/(^-^)z^2_ 



2.5. Transition to Nonrelativistic Evolution 

The transition to nonrelativistic evolution has been discussed for the ISM, constant density 
case by Wijers et al. (1997) and Waxman, Kulkarni, & Frail (1998). For an ambient density 
p = Ar~'^, a straightforward estimate of the transition radius and time is when the shock has 
swept up a mass equivalent to the rest mass of the explosion: 

and 

tNR = ^ = 1.9(1 + z)^yr. (29) 

From equation (^), the relativistic blast wave shock Lorentz factor at this time is F = 1.38 (high 
frequency observation) or F = 1.17 (low frequency observation). 

The nonrelativistic blast wave in a p = Ar^"^ medium is particularly simple for adiabatic 
index ja = 5/3 (e.g., Sedov 1959). The shock radius is given by 

and the postshock profiles are p = ps{r/rs), v = Vgs{r/rs), and p = Ps{r/rs)^, where the subscript s 
refers to the value at the shock front and Vgs is the gas velocity at the shock front. The blast wave 
expansion can be described by = ^o(i?/j4)^/^t^/^, where = 0.78 if 7^ = 5/3 (nonrelativistic 
fluid) and = 0.64 if 7^ = 4/3 (relativistic fluid). The ja = 5/3 evolution yields r^/c = 1.2 at 
the time t-^R. The relativistic and nonrelativistic expansion approximately cross at ^nR) showing 
that this is a reasonable estimate of the transition time. Equation (|29| ) shows that the transition 
to nonrelativistic evolution in a wind is late for typical parameters; the expected value of tNR is 
longer than the typical duration of observations of GRB afterglows. 

In considering the evolution during the nonrelativistic regime, the treatment is similar to that 
in the relativistic regime provided 7m^l- We have R oc t"^^^, shock velocity Vgh oc piv'^f^ oc 
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B cc t ^, and 7^ oc oc t so that oc cx t '''/^ and F^^ cx A^i? cx fit ^ oc t 
where pi is the preshock density and N is the total number of radiating particles. Then, 

(X F^^u^n^'^ oc (z^ < z/„) and F^ oc F^^z^^'^^^^ oc t(5-7p)/6 > example, p = 2, 

yields F^ oc t"^/^ in the nonrelativistic case for u > z^m, but Fy oc in the relativistic case. As 
in interstellar interaction, the light curve steepens after the transition to nonrelativistic flow. 

2.6. Jets 

In most models for GRBs, the burst energy is initially deposited in a relativistic jet. As long 
as the jet is highly relativistic, the observed features should be reproduced by spherical models, 
but as the shocked jet slows there are effects on the observed light curve (as first discussed by 
Rhoads 1997). If 60 is the angular width of the jet, the edge of the jet becomes visible when 
7 K, (Panaitescu & Meszaros 1999). The other important effect is that the slowed jet is able to 
expand laterally. Rhoads (1997, 1999) estimates that this will occur when 7 ~ 0~^/\/3, but Sari, 
Piran, & Halpern (1999) argue that the transition occurs when 7 ~ The issue is the speed 
of sideways expansion and 2-dimensional numerical simulations may be needed to provide reliable 
results. The important point is that during the spreading phase, there is exponential slowing down 
of the forward shock front (Rhoads 1997). Rhoads (1999) has discussed this phenomenon in detail 
for interstellar interaction, but suggests that it also occurs for expansion in a wind. 

Because jet effects are expected when 7 6*"^, equation (^) can be used to find the time 

iM = 2(i^)(^)'i?52^-^day (31) 

when jet effects become important. This expression is appropriate for high frequencies; at low 
frequencies (radio), the coefficient becomes 1 day. A steepening of the afterglow light curve is 
expected at tjet- Once lateral spreading of the jet becomes important, the rapid slowing implies 
that R becomes essentially constant with time and this defines the hydrodynamic evolution (Sari 
et al. 1999). The general nature of the slow expansion is shown by the conservation of energy, 
E ~ c^/3ii?^7^r2, where pi is the preshock density and O is the solid angle of the jet. During the 
jet spreading phase $7 oc 7^ "2 (Rhoads 1997, 1999; Sari et al. 1999) and E « cVl^^ so that R 
must be a weaker function of 7 than any power law, independent of the density profile.^ Thus, 
both ISM and wind interactors should show the same emission properties during the spreading 
phase. We expect F^ oc t^z/^ {v < ua), t'^l'^v^l'^ {va<i' < i^m), t^Pj/^^P"'^^/'^ {u^ < i' < Vc), and 
(j,^ < J,) (Sari et al. 1999). 



^Re'em Sari pointed out this argument. 
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3. PROMPT, REVERSE SHOCK EMISSION 

An initially relativistically hot fireball associated with a GRB cools as it expands, forming 
a cold shell of material coasting at an ultrarelativistic speed when most of the internal energy is 
converted into the bulk kinetic energy. When the coasting shell runs into an ambient medium, the 
bulk kinetic energy is gradually released back into the internal energy, which gives rise to prompt 
emission. Prompt emission has been discussed by Sari & Piran ( 1999a, b) and Meszaros &: Rees 
(1999) in the context of GRB 990123, both assuming a constant ambient density suitable for the 
interstellar medium. In § 3.1, we extend their discussion of prompt emission to the case of a stellar 
wind as the ambient medium. A comparison with the ISM case is presented in § 3.2. 



3.1. Wind Interaction Case 

The structure of the coasting shell is largely unknown, and is determined by the way that 
mass and energy are injected from the central engine. It is likely to be inhomogeneous, since 
bursts of gamma rays are thought to come from one part of the shell running into another. We 
shall, however, ignore such inhomogeneities and adopt uniform distributions for both the shell 
density and Lorentz factor (denoted by 7sh) for simplicity. 

As is well known, the shell interacts with the ambient medium via two shocks: a forward 
shock and a reverse shock. The forward shock runs forward into the ambient medium, whereas 
the reverse shock sweeps up the shell material. The shocked ambient and shell materials are in 
pressure balance and are separated by a contact continuity. Since the thickness of the shocked 
region is expected to be much smaller than its radius, the whole interaction region can be treated 
as planar, as done by Katz (1994) and Sari & Piran (1995). Wc will follow these authors in 
determining the shock properties in a planar geometry, assuming in addition that only a small 
fraction of the internal energy of the shocked materials is radiated away by electrons. 

There are four regions of distinct properties: the unshocked ambient medium (denoted by 
"1"), the shocked ambient medium (denoted by "2"), the shocked shell material (denoted by "3") 
and the unshocked coasting shell (denoted by "4"). We assume that the shocked materials in 
regions 2 and 3 are uniform and move together, and thus share a common bulk Lorentz factor. 
Prom relativistic shock jump conditions, we find that the Lorentz factor is 

712 = (32) 

(measured relative to the nearly static ambient medium, region 1) where ^ = Pa/ Pi is the ratio of 
proper mass densities in the unshocked shell (region 4) and the ambient medium (region 1). The 
proper mass density in the ambient medium is given by 
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where Mw and V„ are the mass loss rate and the velocity of the ambient wind, and R the spherical 
radius. The proper density of the unshocked shell is given by 

where Mgh and A are the proper mass and the width (measured in the frame at rest with respect 
to the origin and the ambient medium and cosmologically local to the burst) of the initial shell. 
The proper mass is related to the total initial kinetic energy Eq of the coasting shell through 

Msh = (35) 

where c is the speed of light. From equations (p3|)-(|35[), we have 



which is independent of radius. Therefore, the shocked materials move at a constant speed 
according to equation (p^). This fact simplifies our discussion of their emission properties. 



The initial shell width A is unknown a priori. It is related to the time Tcr in the frame at rest 
with respect to the origin and cosmologically local to the burst for the reverse shock to cross the 
entire shell through 

= (37) 
c 

The radiation emitted at this time by the shocked shell material along the line of sight is received 
by an observer on Earth at different times, ranging from (1 + z)A/{2c) (for radiation emitted near 
the contact discontinuity) to (1 + z)A/c (for radiation emitted immediately behind the reverse 
shock). The timescale, A/c, therefore characterizes the duration of prompt emission (see Sari & 
Piran 1999b). We assume that the duration is comparable to that of the GRB itself, which is 
typically lO's of seconds for the sources with observed afterglows. That is, the shell needs to have 
a width of tens of light seconds or more. We shall therefore scale A by 10 light seconds, and 
denote the scaled width by Aiq. Scaling other quantities with their typical values, we have 



? = 5-9 , / (38) 



-£^52 
AA1073 

where 73 is the unshocked shell Lorentz factor 7sh divided by 10^ (not to be mistaken as the 
Lorentz factor of region 3). In deriving equation (|37|), we have assumed that ^ <C 73^1 which is 
true for typical parameters. Substituting equation dsl) into equation (|3^), we finally have 



1/4 
52 



^^^ = ^^7W7T74' (39) 

which is much greater than unity unless Aiqj the most uncertain parameter in the above 
expression, is unreasonably large. 
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The shocked ambient medium and the shocked shell material have not only the same bulk 
Lorentz factor, but also the same internal energy density because of pressure balance across the 
contact discontinuity. Relativistic shock jump conditions yield 

63 = 62 = 27sh/9lc2e^/2, (40) 

where the mass density of the ambient wind is given by equation (p^), which can be rewritten into 



with the time T measured in the frame at rest with respect to the origin and cosmologically local 
to the burst. As usual, we assume a (small) fraction eg of the internal energy of the shocked 
matter goes into radiating electrons with a power-law energy distribution. Adopting a constant 
power-law index of p, the minimum electron Lorentz factors in the two shocked regions are 



and 

where nip and nie are the mass of protons and electrons, respectively, and X is the fractional 
abundance of hydrogen, which is close to zero for hydrogen depleted Wolf-Rayet winds. In 
addition, we assume that a constant fraction cb of the internal energy goes into the magnetic field 
in both the shocked ambient medium and the shocked shell material, which yields a total field 
strength of 

„ sl/2 /46bW^^w^V^' 1 .... 

Stot = (Svreijea) / =1 I -, (44) 

in regions 2 and 3. 

To determine the shock emission properties, we adopt the formalism of Sari, Piran & Narayan 
(1998), including cosmological corrections. Let us first consider the reverse shock (region 3). 
Ignoring the synchrotron self- absorption for the moment, there are two characteristic frequencies 
that we need to determine: the "typical" frequency f^s corresponding to the minimum Lorentz 
factor and the cooling frequency Vc- In the observer's frame, we find 

..3 = (^) = 9.8 X 10- ('j^y (^y f^r (^Y' Hz, 

(45) 

where e is the charge of electron and the factor (3p — 6)/(p — 1) equals unity for p = 2.5. We have 
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related the observer's time t to the time T in the frame at rest with respect to the origin through |^ 

l + z T 



2 1I2 



(46) 

The "typical" frequency z^ms shown in equation ( ^Sj ) is to be compared with the cooling frequency 
Uc corresponding to the cooling Lorentz factor given by equation (|To| ) with 7 = 712. The cooling 
frequency observed on Earth is then 



l + zj 2-KmeC ' \l + zj \0.lj 42A^/2 

Clearly, the cooling frequency is much lower than the "typical" frequency for reasonable 
parameters, indicating that most radiating electrons cool quickly down to the cooling Lorentz 
factor 7c. In other words, the shocked shell material is in the fast cooling regime of Sari, Piran & 
Narayan (1998). The emitted flux density therefore peaks at instead of fmS) with a peak value 
given approximately by 

iVe,3(l + ^)P.,max . . 

where the number of radiating electrons N^^^ in the shocked shell region increases linearly with 
time as 

^"'^ (l + z)( 2 ) 7shmpcA ■ 
The peak spectral power -Pj^^max is given roughly by 

Substituting equations (15), (|49|), and ( pO| ) into equation (^8|), we arrive at an approximate peak 
flux at the cooling frequency of 

^ 'V 2 J \i + z - .yTT^. ) \o.iJ 73A10 ^ ' 

which is independent of time. 

With the two characteristic frequencies (i^ms and z/c) and the peak flux (Fj^^^s) thus determined, 
we can now obtain the flux at any given frequency. We are particularly interested in the optical 
prompt emission, say in the R-band at z^^ = 4.5 x 10^^ Hz, a frequency well above the cooling 



"^Strictly speaking, the radiation emitted at time T from the shocked material along the Hne of sight between 
the forward and reverse shocks arrives at the observer at different times. The observed times are (1 + z)T/-fi2, 
(1 + z)r/(27i2), and (1 + 2)7/(4712) for the radiation emitted at time T from the reverse shock front, the contact 
continuity, and the forward shock front, respectively. We pick t = (1 + 2)7/(2712) as a compromise between the 
reverse shock and the forward shock, since the emission from both regions will be computed. 
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frequency Vc but below the typical frequency VrnZ for typical parameters. In this spectral regime, 
we have (cf. Sari et al. 1998) 




,3 



Hc.xif^^^) fer"^l^">Jy. (52) 



1 + Z-.JTT-Z) Ko.il 

Note that in this case the optical flux increases with the square root of the observer's time, and a 
maximum is reached when the reverse shock crosses the inner edge of the freely coasting shell at 
the time 

icr = loi^Aio s. (53) 

The corresponding maximum flux is 



Interestingly, the peak flux is not sensitive to the magnetic energy fraction e^, one of the most 
uncertain parameters. Note that a flux density of 46 mJy in the R-band corresponds to a 
magnitude of 12. Therefore, it is difficult to create a 9*^-magnitude optical flash that lasts for 
several tens of seconds, as observed in GRB 990123 (Akerlof et al. 1999), in the reverse shock if 
the explosion occurs in a typical Wolf-Rayet wind, unless the parameters are far from typical, e.g., 
a much higher explosion energy than 10^^ ergs (see § 4.3 for a discussion). 

For completeness, we now discuss briefly the prompt, optical emission from the forward shock. 
The ambient medium shocked by the forward shock has the same cooling frequency Vc as the shell 
material shocked by the reverse shock, but a substantially higher typical frequency of 



1 A 7i27m2e^- 



tot 



1 + z / InnieC 



fi\ 2 / 1 \2/^ x 1/2 pl/2 

.19 / - O \ / ^ \ / ee2 \ / \ ' 



p-lj \l + Xj \0.lj \0.lj a1/2 



5.7X10- R— : — - ^ Hz, (55) 



t 



where equation (42) is used to eliminate 7m2- For typical parameters, we have Vc < i^r < '^m2, and 
the flux in the R-band is given by 

where the number of radiating electrons A'^g 2 in the shocked ambient medium swept up by the 
forward shock increases linearly with time as 

^,,= fi±^V^-)^|£2ki. (57) 



l + zj rUpV^ 
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Compared with the R-band flux for the reverse shock given in equation (52), the flux in the 
forward shock is down by a factor of 

^ = ^ = J_= 1 A'J'Ali% .... 
F,^,3 iVe,3 2^1/2 4.8 Ell"" ■ ^ ^ 

Therefore, the optical flash is typically dominated by the reverse shock, as in the previously 
studied case of constant-density ambient medium. 

After the coasting shell is completely shocked, the reverse shock front disappears, and there 
is no more shell kinetic energy left to drive the forward shock. Instead of maintaining a constant 
Lorentz factor, the shocked region slows down with time as more ambient medium is swept up. 
The forward shock front begins to evolve as described in § 2. For the material in the reverse shock 
region, we still have Vc < t^m3 at the time that the reverse shock disappears. The electrons in this 
region thus rapidly cool and do not contribute to the emission, which is dominated by the forward 
shock region. 

Synchrotron self-absorption may reduce our estimate of the optical prompt emission. A simple 
way to gauge this effect is to estimate the maximal flux emitted by the shocked shell material as a 
black body 

/ R \ ^ ( 2v^\ 

F,,bb « (1 + ^)'vr j f— jA^Teff, (59) 
(Sari & Piran 1999b) with the observed size given roughly by 

^ 27i2ci, (60) 

and the effective temperature by 

/cTeff 7i27i/"T-eC^/3, (61) 

where is the Lorentz factor of the electrons that radiate at the frequency and is given by 
[27r(l -|- 2;)meCz^/(7i2ei?tot)]^''^- Substituting equations (|60|) and ( |6lD into equation (5£), we have 

l±if { V /"^V"' I'ls-V"* 



which increases rapidly with time as t^/^. At the particular time tcr given in equation (53) when 
the reverse shock crosses the inner edge of the coasting shell, the black body flux is roughly 



^.,.«„^^J (03) 

which is roughly an order of magnitude above the (maximum) prompt optical flux at the same 
time, given by equation (|5^), for typical parameters. Therefore, it appears that synchrotron 
self absorption will not affect our estimate of optical prompt emission from the reverse shock 
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significantly, unless 73 and/or A 10 are unusually small, and/or A-^ and/or £'52 are exceptionally 
large. Finally, we note that the prompt emission estimated here may also be lowered by inverse 
Compton scattering, which cools the radiating electrons in addition to synchrotron emission (Sari 
& Piran 1999b) 



3.2. Comparison with the ISM Case 

Prompt emission from gamma ray bursts in a constant density interstellar medium has been 
investigated in detail by Sari & Piran (1999a, b) and Meszaros &: Rees (1999). For the ease of 
comparison, we shall compute several key quantities of the reverse shock, using the same notations 
and under the same assumptions as the wind case. We limit our discussion to the case of a 
relativistic reverse shock, which demands that the relative Lorentz factor between the shocked and 
unshocked shell materials 

be greater than unity. We assume that the shell 73 is large enough (typically of order unity) so 
that the reverse shock becomes relativistic a few seconds after the explosion. The quantity n in 
the above equation is the number density of the ambient medium in units of cm~^. 

The Lorentz factor for the forward shock is much higher. We find 

712 = 320 , (65) 

which decreases gradually with time. To account for the deceleration, we have adopted a relation 
between the observer's time t and the time T in the rest frame of the origin of t = (1 + z)T/(47^2) 
instead of equation (|4^). Assuming a purely hydrogenic interstellar medium (with X = 1), we 
obtain a typical frequency in the reverse shock of 



-^..^:o»(|^)'X^)(g)^(S)-'%l«- H. 



and a cooling frequency of 



■ N-3/2 A 1/2 

l.exlO^M^) Hz. (67) 



0.1; Eli'nt 

For typical parameters, we therefore have fms < '^r < t'c, and the flux density in the R-band is 
given approximately by 

p f'^msV'-''^/^ 2 y/y 2-V2 y 
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where we have used equations (^)-(^9|) to derive the peak flux density, Fi,^^3,at the typical 
frequency.^ The factor involving the ratio 



Vr Vi + z/ V p- 1 / Vo.i/ Vo.i 

in equation ( |68|) becomes unity for the canonical value of p = 2.5. Note that the optical flux 
density increases with time as t^^^, as in the wind case. For typical parameters, it appears that 
the optical flux in the ISM case is higher than that in the wind case (by a factor of about 6), 
and the emission peaks at a magnitude of about 10 when the reverse shock crosses the inner edge 
of the freely coasting shell. Of course, other choices of parameters can change these numbers 
substantially. For example, if one adopts a shell Lorentz factor 7sh of 300 instead of 10^, then 
the optical fluxes in the wind and ISM cases would become nearly identical, both peaking at a 
magnitude of about 10.7. Note that the prompt optical fluxes in the wind and ISM cases have 
the same dependence on the energy £'52, the shell width Aio, and time t. They have opposite 
dependences on the shell Lorentz factor 73, the ambient density, and the magnetic energy fraction 
€5. Furthermore, whereas the optical flux in the ISM case increases fairly rapidly with the electron 
energy fraction (oc ej if p = 2.5), that in the wind case is independent of eg. 

The X-ray emission from the reverse shock of the wind case is expected to be higher than that 
in the ISM case. For definitiveness, we consider the emission at 2 keV, with = 4.8 x 10^''Hz. 
In the wind case (with X = 0), we have for typical parameters f c < z^^^ < i^ms, and an X-ray flux 
density of 

which increases with time as t^^^. For typical parameters, it has a peak value of 1.4 mJy when 
the reverse shock crosses the inner edge of the freely coasting shell. It is interesting to note that 
the peak value is comparable to the flux density observed in GRB 970228 and GRB 970508, two 
of the GRBs to be discussed in detail in § 4, on the timescale of tens of seconds, which is roughly 
10 keV cm~'^ s~^ (or about 3 mJy at 2 keV; see Figures 1 and 2 of Frontera et al. 1999). The 
prompt X-ray emission from the forward shock should be weaker, as in the case of prompt optical 
emission. Therefore, the prompt X-ray emission in the wind case does not appear to violate any 
observational constraints. Indeed, for some GRBs, it may be able to account for at least part of 
the early X-ray emission detected by the Wide Field Camera of the BeppoSAX. We shall explore 
this possibility elsewhere. 



^We note here that the peak flux would be lower by a factor of about 3 if the formalism of Wijers & Galama 
(1999) is adopted. 
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In the ISM case, we have for typical parameters VmS < I'c < i^x instead, and an X-ray flux 
density of (assuming p = 2.5) 

which is independent of time and is a factor of about 5 lower than the peak flux in the wind case 
for typical parameters. The contrast in prompt X-ray flux between the wind and ISM cases would 
be higher if a smaller shell Lorentz factor is adopted, keeping other parameters at their typical 
values. 



4. OBSERVED SOURCES 
4.1. GRB 970228 

A recent development is the finding of evidence for a supernova in GRB 970228 (Reichart 1999; 
Galama et al. 1999b); the late spectrum is especially convincing. At first sight, this appears to be 
in confiict with the apparently interstellar nature of the afterglow: Fruchter et al. (1999a) found 
time evolution with a = —1.10 it 0.05 at optical wavelengths and the overall evolution appears 
to be compatible with expansion in a constant density medium (Wijers, Rees, & Meszaros 1997). 
However, if the late observations are attributed to a supernova, the decline of the nonthermal 
optical afterglow steepens and a = —1.58 ± 0.28 (Reichart 1999) or a = (Galama et al. 

1999b). With a plausible amount of extinction, Reichart (1999) finds /? = —0.61 it 0.32 at optical 
wavelengths. The optical-to-X-ray spectral index is better constrained; Galama et al. (1999b) find 
PoX = —0.780 lb 0.022 at early times. Cooling evolution {vc below optical wavelengths) is unlikely 
because /3 = —0.78 would imply p = 1.56. With /3 = —0.78, an adiabatic blast wave in a constant 
density medium implies a = —1.17, but a blast wave in a wind implies a = —1.67 (Meszaros et al. 
1998; CL) in good agreement with the observed decline. GRB 970228 can thus be plausibly added 
to the wind interaction category with p ~ 2.6. 

The evolution of the X-ray afterglow shows a = — LSS^Q^g, /3 = —0.96 it 0.19 up to day 4 at 
2-10 keV (Costa et al. 1997) and a = -1.50lo:2| up to day 10 at 0.1-2.4 keV (Frontera et al. 
1998). There is some evidence for a flatter rate of decline at X-ray frequencies compared to the 



optical, which would suggest that I'c is below the X-ray regime. Equation (25) shows that the 
X-ray emission should be in the cooling regime during the observed time period for reasonable 
values of the parameters. The spectrum may be somewhat flatter than /3 = —0.78 in the optical, 
steepening to /? < — 1 in the X-ray regime. 

GRB 970228 was monitored at radio wavelengths, but was not detected at limits between 
10 ^Jy and 1 mJy over the first year (Frail et al. 1998a). In an interstellar interaction model, 
the peak flux, -Fjy^max; should remain constant and move to lower frequency with time. Wijers et 
al. (1997) take -Fi/,max ~ 5 mJy based on the X-ray flux near the end of the initial 7-ray burst. 
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The radio emission clearly did not reach this peak flux. Frail et al. (1998a) explain this result 

by claiming that -F^,max for the afterglow could not be estimated from the initial burst and that 
-fj/jnax ~ 20 — 40 yuJy. This requires a substantial gap in the X-ray evolution between the initial 
burst and the afterglow, and that the first optical observations have occurred shortly after Vm 
passed through optical wavelengths. In a wind model, Fjy^max oc so there is a drop in Fj^^max 
from the time of the early X-ray observations {t ^ 50 s) to the time of the radio observations. The 
radio limits can be accomodated even with an initially high X-ray flux in the afterglow. 



4.2. GRB 970508 

The optical afterglow of GRB 970508 followed power law evolution from day 2 to day ^100 
with a = 1.141 lb 0.014 (Galama et al. 1998a). Galama et al. (1998b) compiled the radio to 
X-ray spectrum of the source on day 12.1 (sec also Wijers &: Galama 1999). They estimated that 
the cooling frequency was at = 1-6 X 10^"^ Hz or a wavelength of 2 fim, i.e. a frequency just 
below optical frequencies. Based on the day 12.1 spectrum, Galama et al. (1998b) estimate that 
p = 2.2, which yields a = —1.15 and /? = —1.1 for u > Uc in both the s = and s = 2 cases. The 
value of p is supported by both the observed spectrum and the rate of decline. However, because 
of the similar optical and X-ray evolution for ISM and wind interaction, radio observations are 
crucial for distinguishing between the s = and s = 2 cases. Extensive radio data exist for GRB 
970508, at 8.46, 4.86, and 1.43 GHz frequencies (Frail et al. 1997; Galama et al. 1998c; Waxman, 
Kulkarni, Sz Frail 1998; Frail, Waxman, & Kulkarni 1999b). Here, we fit these radio data using 
the thin shell model of Li & Chevalier (1999) to determine the burst parameters. Our model 
treats synchrotron emission from a (trans-)relativistic blastwave propagating in an medium. 
Synchrotron self-absorption and relativistic effects are included, but not cooling. As usual, we 
take a power-law distribution of electron energy, and adopt an energy index oi p = 2.2. 

Since the observed radio frequencies are all below the cooling frequency fc in the time interval 
of interest, fitting the radio data alone fixes only three out of four parameters, as discussed in Li 
& Chevalier (1999). The additional constraint comes from the observed fiux in the R-band, which 
is affected by cooling. Together, we find the best model parameters to be: eg = 0.2, = 0.1, 
E52 = 0.3, and = 0.3 (which is in the expected range of a Wolf-Rayet star). The model fit is 
displayed in Figs. 2 and 3; the radio data, from the VLA (Frail et al. 1999b), and the model 
are the same in the two figures and are shown in two different ways for clarity. The general 
evolution shown in the numerical model is along the lines of model D in § 2.4, but it has gradual 
transitions between the evolutionary phases. This is due both to the smoothness of the transitions 
in the radio spectrum and to the time lag effects which result in different parts of the shell being 
observed at different evolutionary phases at the same observer's time. 

The initial rise of the model flux at 1.43 GHz has Fi, oa t, as expected in the self-absorbed 
regime. The model flux on day 6.2 does fall below the 1.43 GHz VLA observation of 100 ± 30 /xJy 
(Frail et al. 1999b). However, monitoring at Westerbork at 1.4 GHz over the first 40 days failed 
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to detect the afterglow (Galama et al. 1998c). The Westerbork observations have lower sensitivity 
than those at the VLA, but a combined map of the 21 observations during this time period yielded 
a flux of 33 lb 40 //Jy. Galama et al. (1998c) take this as evidence that the source is well into the 
self-absorbed regime at early times, as it is in our model. 

Both wind and interstellar relativistic spherical models predict that z/^ oc t^^/^ and the 
data over the range 1.4 — 86 GHz are in approximate accord with this expectation (Frail et al. 
1999b). However, the models differ in that Fy^ decreases with time in the wind case, but not in 
the interstellar case. Fig. 3 shows that the observed decrease of Fy^ is well reproduced by the 
wind model. Additional data are available at 15 and 86 GHz, and Frail et al. (1999b) find that 

^ ly^^^^'^'^'^ : close to the expectation of F^^^^ oc u^^^ in the wind model. 

After about day 100, the 4.86 and 8.46 GHz fluxes tend towards approximate power 
law dechnes (Fig. 2; Frail et al. 1999b). Frail et al. (1999b) find that for t > 90 d, 
ct8.46 GHz = —1.3 lb 0.1 and 04.86 GHz = —1.1 ± 0.1 and that the data at the two frequencies can 
be combined for t > 110 d to yield a = -1.14 ib 0.06 and P = -0.50 ib 0.06 (where F^ oc i°i/^). 
They suggest that "/5 undergoes an abrupt drop from positive to negative values" at t ^ 100 d. 
Our model does not produce an abrupt drop in /?, but we believe that the data are consistent with 
a gradual change (see Fig. 3). We attribute the fact that a4.86 GHz is deduced to be larger than 
^8.46 GHz to the later transition from flat evolution at the lower frequency. Fig. 3 shows that the 
model light curve behavior does go to F^ oc as expected for relativistic, spherical, adiabatic 
expansion in a wind with p = 2.2, but that it eventually becomes steeper. The steepening is due to 
a transition to nonrelativistic expansion, which is included in our numerical model (Li & Chevalier 
1999). At t = 400 d, equation (^) yields 7 = 1.2 for the low frequency case. 

The determination of Uc was carried out as follows. Approximately, the observed late time 
flux decline at 8.46 GHz follows 

/ t \-^-^ 

FsA6 GHz ~ 500 — — fi3y, (72) 



, 100 days 

which implies a flux in the R-band (with i^^ = 4.5 x lO^'^Hz) of 



'8.46 X ion °' / t ^-1-4 



F = Fg 46 GHz ~ 0.73 uJy, (73) 

^ y4.5xl0i4y S.46 Lxiiz V 100 days/ ^ ^ ^ 

if the R-band frequency is below the cooling frequency Vc- The cooling frequency is estimated in 
equation (14) which, for the inferred parameters, becomes 



z., « 1.4 X lOi^tJg Hz. (74) 

Setting Uc = i^j^, we obtain a cooling time = 1000 days for the R-band emission. After t^, 
> Uc, and the R-band flux is given by equation ([7^). Before tj^, we have < Uc, and the 
R-band flux is given instead by 
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This expected flux is shown in panel d of Fig. 1. Our model fits are reasonable overall in all four 

frequencies, lending support to the wind interaction scenario for this GRB. X-ray observations are 
expected to be compatible with our model because it is tlic same as an ISM model for v > V(.. 
The evolution between 6 hr and 6 days can be approximately fitted by Fy oc t~^'^^^'^ (Piro et al. 
1998) and the optical/X-ray spectral index is compatible with (iox = — p/2 = —1.1 (Galama et al. 
1998b). 

Galama et al. (1998b) end their paper advocating an s = blast wave model for GRB 970508 
by describing three deficiencies of the model: (1) -Fi/,max should be constant but is observed to 
decrease with time; (2) va is predicted to be time-independent, but is observed to decrease with 
time, and the rise of the radio fluxes is slower than expected; and (3) the decay after maximum at 
mm wavelengths is perhaps somewhat faster than expected. All of these problems are addressed 
by the wind model, which provides approximate quantitative agreement with the observations. 
Waxman et al. (1998) also noted the problems with a relativistic, spherical, s = model, and 
suggested jet effects and a transition to nonrelativistic expansion as possible solutions. Frail et 
al. (1999b) developed these suggestions in more detail and proposed three phases of evolution 
in a uniform density medium: 1) 2 < t < 25 d, relativistic expansion in which the evolution 
appears spherical (although it is a jet) because relativistic effects allow only a part of the flow 
to be observed; 2) 25 < i < 100 d, jet spreading causes Fy^ and the radio fluxes to drop below 
what would be obtained in an extension of 1); and 3) 100 <t < 450 d, nonrelativistic, spherical 
expansion after jet spreading is complete. Frail et al. (1999b) show that the emission expected 
in phase 3 is consistent with the radio data, but the complete model has yet to be calculated. 
Relativistic time lag effects are likely to be significant at transition times and a possible concern is 
that the optical light curve does not clearly show evidence for the transition from phase 1 to phase 
2. 

Although the radio data appear to provide support for the wind model, the wind interpretation 
does give rise to possible problems regarding synchrotron cooling. In the wind model, Vc increases 
with time; at optical wavelengths, a transition is expected from cooling evolution to adiabatic 
evolution. The opposite is true in s = models. Galama et al. (1998b) cite evidence that Vc 
evolves as expected in an s = model. They identify an observed optical spectral transition 
between 1.0 and 1.8 days from (3 = —0.54 ib 0.14 to /3 = —1.12 it 0.04 with the break frequency Uc 
passing through the Rc band. The wind model would predict cooling in the Rc band at this early 
time and no transition. However, in the time before day 1.5, the observed light curve deviates 
strongly from standard afterglow models (e.g., Galama et al. 1998a; Fruchter et al. 1999b) so 
that the early spectral index cannot be used as a constraint on models for the later power law 
evolution. Galama ct al. (1998b) also cite moderately flat spectra between the K band (Chary 
et al. 1998) and Rc band on days 4.3 and 7.3 as supporting adiabatic evolution in the infrared 
at that time. This is in conflict with the wind model, but we believe that the overall weight of 
evidence supports the wind model. In fact. Chary et al. (1998) find that their data are consistent 
with a t~^ '^ decline and state that their data "agree reasonably well" with a spectral index about 
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-1, as found in the optical (and expected in our model). In the wind model, the optical spectrum 
should flatten and the light curve should steepen at late times. 

Our model fails to account for the optical evolution before day 2 (e.g., Galama et al. 1998a; 
Pruchter et al. 1999b). At early times, the blast wave is expected to be in the fast cooling regime 
{i/m > t'c); this occurs at 500 s in the s = model (Galama et al. 1998b), but at ~ 1 day in the 
wind model. Although the sharp rise in the light curve may be related to this transition, it is not 
possible to account for the rise in a straightforward way. In addition, Piro et al. (1999) have found 
possible evidence for redshifted iron line emission in the X-ray afterglow at an age ~ 1 day. In 
our model, the preshock density at that age is ~ 1 x 10^^^ g cm~'^, which is orders of magnitude 
smaller than the density required for the line feature (Piro et al. 1999). As noted by Piro et al. 
(1999), ordinary stellar mass loss cannot account for the feature that they tentatively observe. 

One expectation of the wind interaction model is that the event be accompanied by a 
supernova. However, Pruchter (1999) has found that a supernova like SN 1998bw added to the 
power law decline of the nonthermal afterglow emission would give an observable bump in the 
light curve at i ~ 20 — 50 days. Such a bump is not seen (Pruchter et al. 1999b) and a supernova 
in GRB 970508 would have to be about 1 magnitude fainter than SN 1998bw (Pruchter 1999). 
Considering possible differences in ejected mass, ^^Ni mass, and explosion energy, we believe that 
some variation in supernova properties is plausible. 

A remarkable property of our model is that it approximately reproduces the radio light curves 
over the entire range of observations from ages of 5 to 400 days. There are reasons why deviations 
from the model might be expected. If the ejecta initially cover a small solid angle, steepening of 
the light curve is expected as the blast wave slows down and there is eventually lateral expansion 
(§ 2.6). Equation (|3T|) shows that 9o^l is required to avoid the steepening; that is, the blast wave 
is spherical or nearly so. Rhoads (1999) and Sari et al. (1999) reached a similar conclusion for 
GRB 970508 based on an interstellar interaction model for the optical emission from the source. 
In addition, the blast wave must remain within the steady stellar wind. Prom equation (^) and 
the blast wave parameters, the observed shock has i? ~ 3 x 10^^ cm. The discussion in CL shows 
that the wind can plausibly extend to this distance. 

4.3. GRB 990123 

The afterglow of GRB 990123 was briefly discussed by CL as a probable case of ISM 
interaction. The observations over the time period 0.01 — 1.5 day (Kulkarni et al. 1999; Galama 
et al. 1999a; Castro-Tirado et al. 1999) can be fitted by an ISM interaction model with p = 2.5 
and optical wavelengths in the adiabatic regime (a = —1.12; f3 = —0.75) and X-ray wavelengths 
in the cooling regime (a = —1.38; P = —1.25). The observed steeper decline in the X-rays vs. the 
R band {ax = -1.44 ± 0.07 vs. or = -1.10 ± 0.03, Kulkarni et al. 1999) is expected for ISM 
interaction, but not for wind interaction. 
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A new feature observed in GRB 990123 was a 9th magnitude optical flash overlapping with 
the GRB (Akerlof et al. 1999). Sari & Piran (1999) and Meszaros k Rees (1999) found that the 
properties of the optical emission over the first ~ 15 minutes could be modeled by synchrotron 
emission from the reverse shock front resulting from interaction with a constant density, interstellar 
medium. In § 3 here, we find it difficult to produce such a strong synchrotron, optical flash with 
standard parameters. If we let z = 1.6, X = 0, and Aio = 4 (so that the emission peaks around 
50 s, as observed), and assume = 73 = 1 and = 0.1, then an explosion energy of 5 x 10^'^ 
ergs is required to produce the 9th magnitude optical flash. Such an energy is nearly two orders of 
magnitude higher than the standard value, but it is not a problem for GRB 990123, the brightest 
7-ray burst with a well localized position. The estimated isotropic 7-ray energy alone for this 
source is ~ 3 x 10^"^ ergs (Kulkarni et al. 1999). The explosion energy could be higher. Therefore, 
the magnitude of the optical flash does not provide a clear discriminator between wind interaction 
and interstellar interaction in this case. 

The observed temporal behavior of the optical flash of GRB 990123 is not compatible with 
the predictions of the wind interaction model. The flash is predicted to rise with time as t^/^ 
(assuming that self-absorption is not significant, which is true for the parameters listed in the 
preceding paragraph), whereas the observed rise is much steeper, close to t^'^ evolution.^ In 
addition, the emission from the reverse shock is short-lived in a wind interaction model and is not 
compatible with the observed flux evolution over the first 15 minutes (Akerlof et al. 1999, Sari 
& Piran 1999a; § 3). The reason for the cut-off of the emission in the wind case is that Vc < ^m, so 
that once the reverse shock front passes through the shell the electrons rapidly cool and there is 
no more emission from the reverse shock. In the ISM case, fc > i^m and long lived emission with a 
power law decline can occur (Sari & Piran 1999a). 

4.4. GRB 990510 

The optical afterglow of GRB 990510 was well-observed in the first 4 days after the burst 
(Stanek et al. 1999; Harrison et al. 1999). The light curve could be well- fitted by an initial 
power law Fi^ cx t"^ followed by a steepening to another power law Fi^ cx . Although there was 
considerable overlap in the data that they used for their analyses, Stanek et al. (1999) found 
ai = -0.76 ± 0.01 and 02 = -2.40 ± 0.02 while Harrison et al. (1999) found qi = -0.82 ± 0.02 
and 02 = —2.18 lb 0.05. The initial flat evolution is strongly suggestive of adiabatic evolution 
(vrn < I' < fc) in a constant density medium (F^^ cx t~'^(P~^)/''). The initial decline rate measured 
by Harrison et al. (1999) is then consistent with p = 2.1 and with the later evolution being due 
to a slowed jet with F^, cx t^^ (Sari et al. 1999). The results of Stanek et al. (1999) suggest a 



* We note here that the steep rise is not reproduced in the simplest ISM case with a uniform freely coasting 
shell and a relativistic reverse shock either. To explain the steep rise, one needs to take into account additional 
complications, such as a nonuniform coasting shell or a transrelativistic reverse shock. 
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somewhat smaller value of p from the early evolution and a somewhat larger value of p from the 

later evolution, but they are close to the expected evolution. The value p = 2.1 implies F^y oc 
which is consistent with the spectrum oc ;/-0-6i±0.i2 nieasured from BVRI photometry (Stanek 
et al. 1999). The hypothesis of jet evolution is supported by radio data which are consistent with 
the expected F^, oc t~^/^ evolution (Harrison et al. 1999). 

Evolution in a constant density medium thus gives a consistent picture for this afterglow. 
Evolution in a wind cannot plausibly account for the early fiat decline, even taking into account 
the possibility of being in the strong cooling regime (§§ 2.2 and 2.3). 

The steepening of the afterglow light curve to a2 = —2.2 after 1.4 days should have facilitated 
the observation of a supernova above the nonthermal emission. Our model would predict the 
absence of such a supernova if it was an interstellar interactor, but a test is difficult because of 
the high redshift of the event, z = 1.619 (Vrccswijk et al. 1999). i/S'T observations by Fruchter 
et al. (1999c) showed that on 17.9 June, 1999 (day 38.5) the afterglow flux was somewhat above 
the extrapolation of the i"^ decline, but the excess counts were a factor 7 below what would be 
expected if a supernova like SN 1998bw were present. 



4.5. Discussion 

A summary of the results from this section is in Table 1. In addition to the sources discussed 
here, we have added GRB 980425 (Kulkarni et al. 1998; Li & Chevalier 1999) and GRBs 980326 
and 980519 (CL). The redshift is given in column 2 for those for which it has been determined (see 
the recent compilation of Ghisellini 1999 for references). The next column gives the afterglow type, 
as determined by the afterglow evolution. We claim that GRBs 970228 and 970508 are probable 
wind interactors, although they have been widely interpreted in terms of interstellar interaction 
in the literature. The reasons for this are varied. In the case of GRB 970228, the recognition of 
supernova emission in the optical light curve yields a steeper nonthermal afterglow decline, in line 
with expectations for wind interaction. In the case of GRB 970508, the cooling frequency Vc is 
somewhat below optical wavelengths so that the high frequency evolution is the same for wind 
and interstellar interaction. However, the predicted evolution is different at radio wavelengths and 
the radio data are compatible with wind interaction. The radio data are crucial for typing this 
burst, and they are suggestive of wind interaction for the afterglow of GRB 980519. 

Although we believe that the radio data for GRB 970508 support a wind interaction model, 
this is a controversial result. Frail et al. (1999b) have proposed a uniform medium interaction 
model, as discussed in § 4.2. In addition, the high frequency, steeply declining afterglows of GRBs 
980326 and 980519 can be interpreted in terms of jet evolution (Sari et al. 1999). If the jet 
model applies, either wind or interstellar interaction could have taken place (see § 2.6). Detailed 
observations over a long time base are needed to clearly discriminate between wind and interstellar 
interaction models. 
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In a paper that appeared after this paper was submitted, Livio & Waxman (1999) addressed 
the problem of distinguishing wind and interstehar interaction afterglows when jets are present. 
They note that once the lateral jet expansion is complete, spherical, sub-relativistic expansion is 
expected (see § 2.5). This should lead to a flattening of the afterglow light curves. If the steep 
decline of some afterglows is due to spherical wind interaction and not to jet interaction, flattening 
of the light curve would not occur unless the blast wave shock front reached the edge of the free 
wind expansion region. 

We have examined the data available on other sources, but have not found more that can 
be typed. The afterglow of GRB 971214 was extensively observed at optical wavelengths, but 
there is probably significant extinction in the host galaxy and radio data are not available. The 
afterglow of GRB 980329 was observed at radio wavelengths; the 8.3 and 4.9 GHz data of Taylor 
et al. (1998) over the first 30 days appear to be in the self-absorbed regime. The data appear 
to be better approximated by Fi, oc t^/^ evolution (interstellar) as opposed to F^y oc t, but there 
are strong scintillation effects. The higher frequency data are difficult to model. Dai &; Lu (1998) 
suggested that GRB 970616 is a wind interactor based on the steep decline with time indicated by 
two X-ray flux measurements; we believe that more detailed observations are needed to make an 
identification. 

Column 4 of Table 1 shows that the estimated electron energy spectral index, p, varies over 
the set of bursts. The probable error in these estimates, ~ 0.1, cannot account for the range. The 
high values of p for GRBs 980326 and 980519 are reduced in jet models for these sources (Sari et 
al. 1999), but a range in p of at least ~ 2.1 — 2.5 is difficult to avoid. The lack of a universal value 
for p calls into question the initial assumption that the shock front accelerates electrons to a power 
law above 7^ that is constant with energy and with time. The possibility of curvature in the 
spectrum and/or time evolution of p should be examined as the data on GRB afterglows improve. 

Column 5 indicates whether there is evidence for a supernova like SN 1998bw in the emission 
from the source. The case is clearest for GRB 980425 and its probable association with SN 1998bw. 
The data on GRB 970228 are compatible with a supernova very much like SN 1998bw (Rcichart 
1999; Galama et al. 1999b). This is also true for GRB 980326, although the observational case is 
not as clear (Bloom et al. 1999b). As discussed in § 4.2, GRB 970508 does not show evidence for a 
supernova like SN 1998bw. Bloom (1999) has reported a similar situation for GRB 980519. These 
facts appear to contradict the expectations of our model. However, in both cases a supernova 
somewhat fainter than SN 1998bw can be accomodated. Supernovae associated with GRBs can 
be expected to have a range of properties; some variation in explosion energy, ejected mass, and 
mass of ^^Ni are all plausible. 

The sixth column lists whether a prompt, optical flash was observed. Such an event has been 
seen only in GRB 990123 and was attributed to synchrotron emission from the reverse shock wave 
(Sari &; Piran 1999a; Meszaros &; Rees 1999). In our picture, the optical synchrotron emission from 
the reverse shock wave of bursts with a wind type afterglow could have a magnitude comparable 
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to, but would die off faster than, that of bursts with an interstellar type afterglow. The decay rate 
of the optical flash in GRB990123 is consistent with an ISM type, not with a wind type. 

The last column of Table 1 lists whether there is evidence for a jet in the afterglow evolution. 
The entry "No" indicates that the jet opening angle is 9o^l- The case of GRB 970508 was 
discussed by Rhoads (1999) and Sari et al. (1999) for ISM interaction and in § 4.2 for wind 
interaction. Li &; Chevalier (1999) found that semi-relativistic, spherically symmetric models 
provided a good description of SN 1998bw/GRB 980425. GRBs 990123 and 990510 show probable 
jet effects on a timescale ~ 1 — 2 days (Kulkarni et al. 1999; Harrison et al. 1999; Sari et al. 1999). 
The cases of GRBs 980326 and 980519 are controversial. If the steep declines of their optical 
afterglows are interpreted as jet effects (Sari et al. 1999) then the jet effects occur at an earlier 
time than for GRBs 990123 and 990510, but if the declines are interpreted as wind interaction 
(CL) then the jet effects occur later if at all. 



5. TWO TYPES OF PROGENITORS 

In CL, we suggested that there are two types of GRB progenitors among the well-observed 
sources: massive stars that have afterglows characteristic of wind interaction and are likely to be 
accompanied by supernovae, and compact binary mergers that have afterglows characteristic of 
constant density, interstellar interaction and are not accompanied by supernovae. Our discussion 
here has increased the wind group by adding GRBs 970228 and 970508. The pattern of association 
with a supernova is strengthened by the finding of strong evidence for a supernova in GRB 970228. 
A new result here is that the prompt optical flash in a wind interaction GRB is expected to rise 
slowly and to disappear abruptly, in contrast to the observed temporal behavior of the optical 
flash of GRB 990123. The contradiction points to an interstellar interaction for this case, which 
supports the interpretation of the afterglow light curve evolution by CL. 

Another possible indicator of the GRB progenitor type is the location in a galaxy (Paczyhski 
1998). For the cases in which imaging with HST is available, the bursts appear to be superposed 
on or near the optical disks of the host galaxies. In the case of GRB 990123, the burst is at about 
5.8 kpc (0."67) from the center of the host galaxy, near the edge of the optical disk (Bloom et al. 
1999a). In the case of GRB 990510, no host galaxy has been detected, even with deep imaging 
(Israel et al. 1999; Pruchter et al. 1999c). However, the production of a typical afterglow for the 
interstellar interaction case probably requires that the burst take place in or near a galaxy so that 
the surrounding density is sufficiently high. For GRB 970508, which we have identified as a wind 
interactor, the burst is within 0."01, or 70 pc of the galaxy center (Fruchter et al. 1999b). Wind 
interaction implies a massive star progenitor and the source may have been in a nuclear starburst 
region. In the case of GRB 970228, the source is about 0."5 from the galaxy center (Sahu et al. 
1997), which is comparable to the interstellar interaction case of GRB 990123. Location does not 
appear to provide a clear discriminator for the two types of bursts considered here. 



- 30 - 



We have argued that there are two types of burst progenitors, so the GRBs associated with 
the two types might be expected to have distinct properties. We have examined the data on the 
GRBs themselves and have found no clear distinction between those that give rise to wind or to 
ISM type afterglows. In both cases, the 7-ray light curves are complex and last for lO's of seconds. 
One implication is that the mechanism giving rise to the initial GRBs does not depend on the 
external medium. This supports the generally favored internal shock model for the GRBs (e.g., 
Piran 1999). 

Another implication is that, for both types, the source must be capable of producing a long 
duration burst. This is expected for a massive star origin because the collapse time of the stellar 
core is on the order of lO's of seconds (MacFadyen & Woosley 1999; Fryer, Woosley, & Hartmann 
1999a). Among the compact binary mergers that have been proposed (Fryer et al. 1999a and 
references therein), the type that may satisfy the duration requirement is the black hole - white 
dwarf merger (Fryer et al. 1999b). In addition, this type of merger is expected to occur in or 
near the host galaxy, again in accord with observations. For the small number of sources we have 
discussed here, the rate of massive star explosions is not very different from the rate of compact 
star explosions. However, the observed interstellar interactors are somewhat more distant and 
luminous than the wind interactors on average, and the occurrence of GRB 980425/SN 1998bw 
indicates the presence of low energy wind interactors. This would imply a higher rate of massive 
star explosions per unit volume than of compact star explosions. However, if jet effects are more 
important for interstellar interactors (as suggested by Table 1), the rate of interstellar interactors 
is increased, so no clear conclusions are possible for the relative rates. 

One of the unexpected results of this study is the possible evidence for a relation between 
jet effects and the afterglow (and thus progenitor) type. Based on a small number of objects, 
the interstellar interactors show evidence for significant jet effects, while the wind interactors 
with massive star progenitors do not. This suggests the intriguing possibility that the GRB 
engine generates a collimated flow, but in the case of a massive star progenitor the flow becomes 
uncollimated upon passing through the star. 

6. CONCLUSIONS 

The results of our work can be summarized as follows: 

1. The afterglow light curves resulting from the interaction of a GRB with a circumstellar 
wind have properties that depend on the observation frequency. At radio frequencies, there are 
two possibilities. Above VAm (equation @|), the light curve is of type D (Fig. 1) and below UAm, 
it is of type E. Detailed integrations over a spherical blast wave show these segments, but without 
sharp transitions. At high frequencies (optical and X-ray), the expected light curve is of type A 
(Fig. 1). The evolution may make a transition from cooling [Fi, oc v' 

-p/2^-(3p-2)/4) adiabatic 

{Fp oc zy~(P~^)/2i-(3p-i)/4 ^ evolution, or may remain in one of these phases during the period of 
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observation depending on the parameters, especially es- The transition occurs earlier at lower 
frequency. For typical parameters, the fast cooling phase (all electrons cool) lasts for ~ 1 — 2 days, 
considerably longer than in the interstellar interaction case. The light curve can be modified from 
the above expressions during this early phase. 

2. During wind interaction, the transition to nonrclativistic evolution occurs at ~ 2 yr for 
typical parameters. This is longer than in the interstellar interaction case because of the low 
density in the outer parts of the wind. 

3. The optical synchrotron emission from the early reverse shock wave has a peak magnitude 
of ^ 12 for standard parameters. It could be comparable to the prompt emission for the interstellar 
interaction case. Radio emission is strongly self-absorbed during this phase. 

4. The recognition of supernova emission in the optical radiation from GRB 970228 (Reichart 
1999; Galama et al. 1999b) implies that the decline of the nonthermal afterglow is faster than 
previously thought and is compatible with adiabatic wind interaction. 

5. The optical spectrum and decline over days 2 — 100 of GRB 970508 suggest cooling 
evolution and are compatible with either interstellar or wind interaction models. However, the 
extensive radio data on this source strongly suggest wind interaction. In our model, the blast wave 
energy is 3 x 10^^ ergs and the wind mass loss rate is 3 x 10~^ Mq yT~^ for a wind velocity of 
1,000 km s-^ 

6. In the case of GRB 990123, the afterglow showed a steeper decline with time in X-rays 
than in optical emission, as expected for interaction with a constant density medium. A prompt 
optical flash was observed from this object, which can be plausibly attributed to synchrotron 
emission from the reverse shock front due to interaction with the interstellar medium (Sari Sz Piran 
1999a,b). We have found that interaction with a wind cannot produce the temporal behavior 

of the observed flash. The case for interstellar interaction appears to be strong. We have also 
identified GRB 990510 as a likely interstellar interactor based on its early time evolution. 

7. Among the 7 sources for which we have suggested identifications, 5 are wind interactors 
and 2 are interstellar interactors, although two of the wind cases can also be interpreted as jets in 
either a wind or the interstellar medium. The interstellar interactors are somewhat more luminous 
and distant, but appear to have more significant jet effects, so no clear conclusions can be drawn 
regarding the relative rates of the two types. However, there does not appear to be a clear 
distinction in the initial 7-ray burst properties of these different types. This provides additional 
evidence for a burst model that does not depend on environment, e.g., the internal shock model 
for bursts. 

The main result of our paper is the presentation of evidence for two types of GRB afterglows 
in different environments: a constant density interstellar medium and the wind of a Wolf-Rayet 
star. The types are not immediately distinguishable because at an age of a few days, the preshock 
wind density is comparable to an interstellar density. At an age of seconds, the preshock density 
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is higher for the wind case and we predict that the wind interactors have prompt optical emission 
that die off faster than the interstehar interactors. In addition, we expect the wind interactors to 
be accompanied by a supernova event. Future observations of GRBs and their afterglows should 
provide a clear test of our model. 

We are especially grateful to the referee, Re'em Sari, for pointing out corrections and for 
constructive comments in his detailed report. Comments and correspondence from Eli Waxman 
and Dale Frail were important for clarifying and correcting our model for GRB 970508, although 
they do not agree with our model. We also thank Andy Fruchter and Pawan Kumar for useful 
correspondence and conversations. Support for this work was provided in part by NASA grant 
NAG5-8232. 
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Table 1. Properties of GRB Afterglows 



Burst 


Redshift 


Afterglow 


Spectral 


Supernova 


Prompt, Bright 


Jet 


GRB 


z 


Type 


Index, "p 


like SN 1998bw 


Optical Flash 




970228 


0.695 


Wind 


2.6 


Yes 






970508 


0.835 


Wind 


2.2 


No 




No 


980326 




Wind 


3.0 


Yes 










Jet 


2.2 


Yes 




Yes 


980425 


0.0085 


Wind 


2.5 


Yes 




No 


980519 




Wind 


3.0 


No 










Jet 


2.2 


No 




Yes 


990123 


1.60 


ISM 


2.5 




Yes 


Yes 


990510 


1.619 


ISM 


2.1 


No 




Yes 
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Fig. 1. — Characteristic light curves for wind interaction models in various frequency ranges. The 
light curves are sorted from high frequency at the top (A) to low frequency at the bottom (E). 
X-ray and optical light curves are typically of type A and radio of type D. 
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Fig. 2. — Wind interaction model for the afterglow of GRB 970508. 
Frail et al. (1997; 1999b). R-band data are taken from Sokolov et al. 
parameters are listed in the text. 



Radio data are taken from 
(1998). The best model fit 
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Fig. 3. — Log-log plot of the same observed data and model fits at the radio frequencies as in 
Figure 2, showing the curvature in the power-law flux decay after about day 100. The data and 
fits for 4.86 GHz and 1.43 GHz are artiflcially lowered by a factor of 10 and 10^, respectively, for 
clarity. 



